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We study the evolution of time-dependent fluctuations and particle production in an expanding dS 
and contracting AdS universe. Using the functional Schrodinger formalism we are able to probe the 
time dependent regime which is out of the reach of the standard approximations like the Bogolyubov 
method. In both cases, the evolution of fluctuations is governed by the harmonic oscillator equation 
with time dependent frequency. In the case of an expanding dS universe we explicitly show that 
the frequency of fluctuations produced at a certain moment diminish in time, while the distribution 
of the created particles quickly approaches the thermal radiation of the dS space. In the case of a 
contracting AdS universe we show that the frequency of fluctuations produced at a certain moment 
grow in time. Nominally, the temperature of radiation diverges as the Big Crunch is approaching, 
however, increasing oscillations of the spectrum make the temperature poorly defined, which is in 
agreement with the fact that AdS space does not have an event horizon which would cause thermal 
radiation. Unlimited growth of fluctuations indicates that an eventual tunneling into AdS vacuum 
would have catastrophic consequences for our universe. 



I. INTRODUCTION 



Space-times with the constant vacuum energy density 
play very important role in modern cosmology. It is be- 
lieved that the primordial inflation was driven by the 
cosmological constant or a scalar field vacuum energy 
density of positive sign. As the recent observational data 
indicate, our late time universe may again be dominated 
by the (positive) vacuum energy density which drives the 
accelerated expansion of the universe. Both primordial 
and late time phases of acceleration can be effectively de- 
scribed by de Sitter (dS) space-time. On the other hand, 
as indicated by string theory and other generic consider- 
ations PJ, it may happen that the true vacuum energy 
density of our universe is negative, which is described by 
anti-de Sitter (AdS) space. In that case, our universe will 
in the future undergo a phase transition from dS to AdS 
space. 

dS space describes an expanding universe. Time de- 
pendent gravitation field usually leads to particle pro- 
duction. Particles that propagate in the background of 
a time dependent metric of an expanding universe get 
excited and radiation is produced. For example, this is 
the mechanism behind particle production in FRW uni- 
verse. In addition to this, dS space contains an event 
horizon (called dS horizon, not to be confused with the 
particle horizon). Using an analogy with the event hori- 
zons in black hole space-times, one could expect that 
the produced radiation of dS space is thermal with the 
constant temperature which is proportional to the value 
of the cosmological constant A. Indeed, such a result 
may be obtained using standard approximations like the 
Bogolyubov method or tunneling formalism. However, 
time-dependent evolution of radiation is beyond the reach 
of these standard approximations. The main goal of this 



paper is to study these time-dependent effects. We will 
use the recently developed functional Schrodinger formal- 
ism to study a scalar field propagating in the back- 
ground of an expanding dS space. We will show that the 
evolution of fluctuations (excitations) of the scalar field 
is governed by the harmonic oscillator equation with the 
time dependent frequency. We will solve the equations of 
motion to find the time-dependent wave function for the 
system. If we set the wave function of the system to be 
in the ground state at a certain moment of time, the par- 
ticle content at some later time will be given by the wave 
function overlap between the initial and final state. This 
wave function overlap will give us the time-dependent 
occupation number of produced particles. In cosmology 
one is very often interested in the power spectrum, i.e. in 
power contained in a mode of a certain frequency, but this 
quantity is very closely related to the occupation number 
of particles in a given mode Q . We will show that when a 
mode of ceratin frequency is created it does not automat- 
ically follow a thermal distribution with the temperature 
proportional to the value of the Hubble parameter (which 
is at late times proportional to the cosmological constant 
A). Modes whose frequencies are much larger or much 
smaller than the Hubble parameter do not follow a ther- 
mal distribution, while those of the order of the Hubble 
parameter do. As the high frequency modes get stretched 
by expansion, and their frequencies become comparable 
to the Hubble parameter, they become thermal. 

We will then apply the same technique to AdS space. 
In the cosmological context, AdS space may be used to 
describe a collapsing universe. Unlike dS space, AdS 
space does not contain an event horizon, so we do not 
expect thermal radiation with the constant temperature. 
However, a time-dependent metric may again lead to par- 
ticle production. The fate of the universe which tunnels 
into the true vacuum which is AdS was studied in [J, u3f . 
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There, it was argued that after a short period of expan- 
sion, the universe starts recollapsing and ends up form- 
ing a black hole. Here, we study the time-evolution of 
the scalar field propagating in the background of AdS 
space and our findings support these arguments. Fluctu- 
ations of the scalar field produce more and more modes 
of high frequency at late times. The best fit temperature 
of the produced particles diverges at the time of the Big 
Crunch. However, towards the Big Crunch, oscillations 
of the spectrum grow in amplitude making the tempera- 
ture less and less defined. 



II. FRIEDMAN-ROBERTSON-WALKER 
SPACE-TIME 

To setup the formalism, we will first consider the stan- 
dard Friedman-Robertson- Walker (FRW) space-time. 
The action for the scalar field propagating in the curved 
space-time background given by the metric tensor g^ v is 



(1) 



Decomposing the (spherically symmetric) scalar field into 
a complete set of real basis function denoted by {f k (r)} 



(2) 



we want to find a complete set of independent cigenmodes 
{b k } (which are linear combinations of the original modes 
{a k }) for which the Hamiltonian is a simple sum of the 
independent modes. The total wavefunction then fac- 
torizes and can be found by solving a time-dependent 
Schrodinger equation of just one variable. 

The metric for the spherically symmetric FRW space- 
time is 



ds 2 = -dt 2 + R{tf 



dr 2 



1 - kr 2 



'-dn- 



(3) 



where 



dD. 2 = d9 2 + sin' 



R(t) is the scale factor, while k 
and Eq. (TT]) we can write 



S = 2ir / dt / drr 



VI - kr 2 



(4) 

1,0,1. Using Eq. © 
1 - kr 2 



(5)' 

Using the mode expansion in Eq. @, we can rewrite 
the action as 



S 



dt 



~p-a(t) k Mkk'Ci{t)k> + ^Y-a(t) k Nkk'a(t)k> 

(6) 



where d(t) — da(t)/dt, and M and N are matrices that 
are independent of R(t) and are given by 



M fcfc , = 4tt / drr 



1 



ffkfk' j 



N kk , = 4nf drr 2 ^l-kr 2 f k f' k ,. (8) 

For our procedure we will use the fact that the matri- 
ces M and N are symmetric, while their exact form will 
not be needed. From the action Eq. ([5]) we can find the 
Hamiltonian, H, and write the Schrodinger equation 



H(a k ,t)ip(a k ,t) = i 



.dip(a k ,t) 



c)t 



(9) 



The wavefunction ip{a k ,t) is actually the wavcfunctional 
since a k (t) is a function itself. Therefore, this procedure 
is the functional Schrodinger formalism. Defining n^, the 
momentum operator conjugate to a(t) k as 



n fc = -i- 



d 



da(t) k 

the Schrodinger equation can be written as 
1 „ „ R(t) 



2R 3 (t) 



-a(t) k N kk ,a{t)v 



(10) 



(11) 

The principal axis transformation guaranties that M and 
N can be simultaneously diagonalized. Therefore, the 
Schrodinger equation Eq. (fTTj) decouples into an infi- 
nite set of decoupled equations. The single eigenmode 
Schrodinger equation is 



1 d 2 K , . w2 



W,t)=iRHt) d -^ (12) 



where m and K denote eigenvalues of M and N, and b 
is the eigenmode (which is a linear combination of the 
original modes a k ). 

Re- writing Eq. (fT2")) in the standard harmonic oscillator 
form (with the time dependent frequency) we obtain 



_L_ cP_ 



where 



and 



-uj 2 (r,)b 2 



K 



tp(b,rj) 



.dip(b, rf) 



u\r,) = -R*{t) = utR*{t) 



rj = 



dt' 



R 3 {t>) 



(13) 



(14) 



(15) 



Obviously, ujq = y ^ is the initial frequency of the mode 
at the time of its creation. The exact solution to Eq. (fT3")) 
is given by 0] 



i/>(b, V ) = e ta M ( ^ 
7rp 



1/4 



cxp 



(16) 
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where = dp/dr] and p is given by the real solution of 
the ordinary (though non-linear) differential equation 



p m + uj 2 (i])p = i 



with initial conditions 
p(0) = 
The phase a is given by 



, ^(o) = o. 



1 r 



dij 



(17) 



(18) 



(19) 



Consider an observer with detectors that are designed 
to register particles of different frequencies for the free 
field $. Such an observer will interpret the wavefunction 
of a given mode b at some later time in terms of simple 
harmonic oscillator states, {<p n }, at the final frequency, 
ui. The initial (t = 0) vacuum state for the modes is the 
simple harmonic oscillator ground state 



<p(b) 



(20) 



The number of quanta in eigenmode b can be evaluated 
by decomposing Eq. (|16[) in terms of the states and eval- 
uating the occupation number of that mode. Writing the 
wavefunction for a given mode in terms of simple har- 
monic oscillator basis at t = is given by 



where 



^(6,t) = 5^c»(t) Vn (6) 



(21) 



(22) 



which is an overlap of a Gaussian with the simple har- 
monic oscillator basis functions. The occupation number 
at eigcnfrcqucncy ui (which is the frequency at the final 
time of the measurement corresponding to tf), is given 
by 



N(t,Ul) = J2 n \Cn\ 



(23) 



The occupation number in the eigenmode b is then 
given by (see Appendix |A"| 



N(t,ui) 



up 

7f 



i 



>jjp i 



En. 

Ulp 



(24) 



This occupation number represents the cumulative num- 
ber of particles produced by the time tf. To extract 
the temperature corresponding to the occupation number 
distribution, we can compare Eq. (|24]l with the Planck 
distribution 



N P (cu) = 



1 



- 1 

where j3 is the inverse temperature. 



(25) 



III. PARTICLE PRODUCTION IN DE SITTER 
SPACE-TIME 

We now consider particle production in the background 
of dS space-time. The metric for dS space in static coor- 
dinates can be written as 



ds 2 



1 



r 



dr 



i - 



r 2 



dr 2 



2 dQ 2 . (26) 



The constant I is called dS radius. Obviously, at r = I 
the metric has an event horizon. This horizon is not to 
be confused with a particle horizon. The particle horizon 
represents the largest comoving distance from which light 
could have reached us by now, while the event horizon 
is the largest comoving distance from which light emit- 
ted now can ever reach the observer at any time in the 
future. Since black hole solutions with an event horizon 
are known to produce thermal radiation, it is expected 
that dS space can produce thermal radiation as well (for 
some counterarguments see 0,0). 

We can also write the metric for dS space in non-static 
coordinates. For this purpose, we set k — in Eq. ([3]). 
The solution to the Einstein's equations with a positive 
cosmological constant X ^ l/l give the expansion factor 
R(t) ~ exp(Ht), with H (the Hubble parameter) given 
late times by 

H = y/X/3. (27) 

Therefore, the time dependent dS metric is given by 

ds 2 = -dt 2 + R 2 (t) (dr 2 +r 2 dn 2 ) , (28) 

with R(i) ~ exp(Ht). From here, we proceed like in 
Section HU Since k = 0, the matrices M and N, Eqs. ([7]) 
and JHI) are now 



M = 4tt / drr 2 hh 



N 



4n / drr 2 f k f k ,. 



(29) 
(30) 



The occupation number is given again by Eq. (|24[) . 

We comment here on variety of frequencies we use 
in this context. A frequency certainly depends of the 
time evolution parameter that an observer is using. The 
Schrodinger equation Eq. (|13|) is written in terms of 77, so 
ui is the mode frequency with respect to rj and not with 
respect to time t (which we denote with u/*'). From 
Eq. (fTS"]) we see that the frequency in t is R(t)~ 3 times 
the frequency in rj. Since TV is expressed in terms of Q, 
which is the frequency at the final time corresponding to 
tf, we have 



-,(*> 



-3Ht f 



(31) 



From Eq. (fl"4|) and Eq. ([3T|) we see that the final fre- 
quency of the harmonic oscillator states is 



u> 



(0 



w e 



-Ht f 



(32) 
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ZZ^/H = 2, 3, 4 

N 




FIG. 1: The occupation number of produced particles, N, 
as a function of dimensionless time, tH, for various values of 
the (final) dimensionless frequency uj^ /H, in dS space. The 
curves are higher for lower values of ui^/H (higher wave- 
lengths). At some finite time, the occupation number reaches 
the plateau and apart for some diminishing fluctuations re- 
mains constant at late times. 

This means that the frequency of the mode whose ini- 
tial frequency was uiq decreases with time, as expected. 
The physical frequency measured by the detector at some 
time tf is Q^>. Expression ([52)1 says that a ceratin fixed 
value of ui^\ at some late time, comes from a mode 
of much higher original frequency loq, which was sub- 
sequently stretched by the expansion. 

In Fig. [1] we plot the occupation number of produced 
particles, N, as a function of dimensionless time, tH , for 
different values of dimensionless frequency u)^> /H. The 
curves are higher for lower values of u>^ / H, which means 
that the occupation number in lower frequencies (higher 
wavelengths) is higher. At some finite time, the occu- 
pation number reaches the plateau and apart for some 
diminishing fluctuations remains constant at late times. 
This can be explained with the help of Eq. (j3"2")l . Modes in 
an infinite range of frequencies are produced. The modes 
of any particular frequency get stretched by the expan- 
sion, but they also get constantly replaced by stretched 
modes whose initial frequency was smaller. 

Note that / H is the frequency of the harmonic os- 
cillator at the final time tf at which the occupation num- 
ber is evaluated. With fixed tf, in order to vary Q^/H 
we need to vary ujq since ui^ = ujQe~ Ht . Thus, different 
values of ui^/H for the fixed t = tf, correspond to dif- 
ferent wo which is the original frequency of the mode at 
the time of its creation. 

In Fig. [2J we numerically evaluated the spectrum, N vs. 
ujW/H, of mode occupation numbers at several values of 
tH . The curves are higher for later times, which means 
that there are more excited modes at later times. In 
agreement with the plot in Fig. [TJ at any fixed time there 
are more lower frequency modes. As expected, at later 
times we have progressively more lower frequency (longer 



tH = 0.75, 1, 1.5 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 ' 

FIG. 2: The occupation number of produced particles, N, as 
a function of the dimensionless frequency Q^' /H, for various 
values of dimensionless time, tH, in dS space. The curves 
are higher (contain more particles) for later times. The plot 
shows the non-thermal features at very low and very high 
frequencies. 

wavelength) modes produced by the expansion. 

Comparing the distribution in Fig. [2] with the ther- 
mal Planckian distribution Eq. (|25[) we can infer that 
that only the modes with frequencies comparable to the 
inverse Hubble radius follow the thermal distribution. 
Those with much lower and much higher frequencies show 
departure from the thermal distribution. The thermal 
distribution in Eq. (|25[) is divergent at low frequencies 
while the distribution in Fig. [5] is not. At high frequen- 
cies, the distribution in Fig.[2]has some non-thermal oscil- 
lating features. Short wavelength modes become thermal 
when the expansion stretches them to the scales compa- 
rable to the Hubble radius. 

In order to fit the temperature that corresponds to the 
particle distribution, we first express N(u>) in terms of 
time t. Then we plot ln(l + 1/N) versus ui^ / H and find 
the slope /3. Finally, we correct for the factor in Eq. (f5Tj) 
which implies 

T = e - 3m fj3- 1 (tf). (33) 

This is then the temperature measured at some moment 
tf, seen by the observer who is using time t as an evolu- 
tion parameter. 

In Fig. [3] we plot ln(l + 1/N) versus Q^' / H for various 
values of tH. Here we see that the curves corresponding 
to tH = 3 (blue) and tH = 4 (red) are almost indistin- 
guishable. Since ln(l + 1/N) versus ui^'/H will give us 
information about the temperature, this means that the 
temperature does not change at late times. 

Let us analyze the plot in Fig. [3] We mentioned that 
only the modes with frequencies comparable to the Hub- 
ble constant (inverse Hubble radius) follow the thermal 
distribution. Therefore, in order to get the temperature 
of the produced radiation, we fit the straight line only 
in the segment of the plot corresponding to the neigh- 
borhood of uW/fl = 1. From Fig. H we see that most 
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tH = 3, 4 

Ln(l + 1/N) 




FIG. 3: Here we plot ln(l + l/N) versus lj^/H for the val- 
ues of tH = 3, 4 for the dS space. The curves for tH = 3 
(black) and tH = 4 (red) are almost indistinguishable which 
indicates that the temperature does not change with time. 
Since only the modes with /H ~ 1 are thermal, we fit the 
temperature only in that segment of the plot and find it to 
be proportional to the Hubble constant. Very high and very 
low frequencies are not thermal. 

of the modes belong to that range (say we sum up the 
number of modes within one order of magnitude, from 
ui^/H = 0.15 to ui^/H = 1.5). The number of modes 
of very high frequency is very low since Nfa^') drops in 
that regime quickly, while the number of modes with very 
low frequency is very low since that regime is a very short 
segment on the N(ui^) curve. Thus fitting the curve in 
the omega ui&'/H ~ 1 segment would give a pretty re- 
alistic result for the temperature of the produced radia- 
tion since most of the produced particles will be in that 
regime. Applying Eq. ([55]). the fit gives for the slope 
/3 1 in units of 1/ff which implies 

T ~ if. (34) 

As expected, the temperature of the produced radiation 
is proportional to the Hubble constant, which is in turn 
proportional to the value of the cosmological constant A. 

The modes with much lower and much higher frequen- 
cies than the Hubble constant do not follow the ther- 
mal distribution. However, at some later time, modes 
that were created with short wavelengths become ther- 
mal when the expansion stretches them to the scales com- 
parable to the Hubble radius. 

IV. PARTICLE PRODUCTION IN ANTI-DE 
SITTER SPACE-TIME 

Now we apply the analogous study to AdS space with 
the negative cosmological constant, i.e. A < 0. For this 
purpose, we set k = — 1 in Eq. ([3]). In that case, the solu- 
tion for the scale factor is given by R(t) — if -1 cos(ifi) 
where if is if = ^/|A|/3. Therefore the time dependent 
AdS metric takes on the same form as dS metric given 



in Eq. ([28]) , apart from the different scale factor R(t). 
For the decaying segment of cos(ifi), this time depen- 
dent solution describes a collapsing universe under the 
influence of the negative constant vacuum energy den- 
sity represented by the negative cosmological constant. 
The matrices M and N are now given by Eq. (|2U)) and 
Eq. (f3"0")l . with the substitution of k = —1 

M kk , = 4tt f drr 2 -L= f k f k ,, (35) 
J Vl + r z 

N kk , = 4tt J drrWl + r 2 f k f' k ,. (36) 

The occupation number is given again by Eq. (|24l) . 
Repeating the procedure from the previous section, we 
plot N versus tH for various values of ui^ /H in Fig. |4] 
The Figure shows that for different values of dimension- 
less frequency Co^/H, the occupation number increases 
in all frequency modes. However, the curves are higher 
for higher values of u)^ /if, which means that the occu- 
pation number in higher frequencies (lower wavelengths) 
increases faster at later times. In contrast with dS space 
where the occupation number at some fixed final fre- 
quency is constant at late times, the occupation number 
in AdS space diverges, and it does this in finite time, 
indicating the big crunch. 



ZZ^/H = 1, 2, 5 

N 
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FIG. 4: The occupation number of produced particles, N, 
as a function of dimensionless time, tH, for various values of 
dimensionless frequency cj^'/H, in AdS space. The curves are 
higher for higher values of w'*'/if. The occupation number 
diverges in finite time, indicating the Big Crunch. 

In Fig.[5j we have numerically evaluated the spectrum, 
N versus u)W /if, for several values of tH. The curves are 
higher for later times, which means that there are more 
excited modes at later times. As expected, at later times 
we have progressively more higher frequency modes since 
the space is collapsing. 

We may now try to fit the temperature that corre- 
sponds to the particle distribution. We first express N(u>) 
in terms of time t and then plot ln(l + l/iV) versus ui^ /if 
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for various values of tH . We could then fit the straight 
line through the curves and find the slope j3 for each of 
them. Eq. (fT5|) tells us that the frequency in t is R{t) 3 
times the frequency in 77, and at time tf, this implies 



-,(') 



H' 



cos 3 (Ht) 



(37) 



Therefore the temperature seen by the observer whose 
parameter of evolution is t should be 



T 



H> 



cos 3 {Ht) 



(38) 



However, following this procedure we find it very dif- 
ficult to precisely fit the temperature. Nominally, from 
Eq. (|38|) , the temperature diverges as tH approaches the 
time of the Big Crunch. However, from Fig. [5] we see 
that as tH increases the oscillations of the curves increase 
making the temperature less and less defined toward the 
Big Crunch. This is in agreement with the fact that AdS 
space (unlike dS space) does not have an event horizon 
which could cause thermal radiation. 



V. CONCLUSION 

In this paper, we study the time-dependent evolution 
of the scalar field propagating in the backgrounds of an 
expanding dS and contracting AdS universe. Using the 
functional Schrodinger formalism we are able to probe 
the time dependent regime which is out of the reach of the 
standard approximations like the Bogolyubov method. 
In both cases, the evolution of fluctuations is governed by 
the harmonic oscillator equation with time dependent fre- 
quency. In general, the harmonic oscillator equation gov- 
erns the evolution of quantum fluctuations and particle 



N 



tH = 1.4, 1.45, 1.5 




FIG. 5: The occupation number of produced particles, N, as a 
function of dimensionless frequency d/'' / H, for various values 
of dimensionless time, tH, in AdS space. The occupation 
number at all frequencies grows as tH increases, however the 
high frequency modes grow at a higher rate. 



Ln(l + 1/N) 



tH = 0.5,tt/4, 1, 1.5 




FIG. 6: Here we plot ln(l + l/N) versus cd'*' /H for various 
values of tH for the AdS space. The curves are lower for later 
times. The slope of the best fit straight line is supposed to 
give the temperature of radiation. However, the plot shows 
that the oscillations become greater as tH increases, making 
the temperature less and less defined toward the Big Crunch. 



production. We solved the equations of motion and found 
the time-dependent wave function for the system. The 
wave function overlap between the initial state (which is 
the vacuum) and the final state at some arbitrary later 
time gave us the time-dependent occupation number of 
produced particles. 

In the case of an expanding dS universe we explicitly 
showed that the frequency of fluctuations produced at a 
certain moment diminish in time, while the distribution 
of the created particles quickly approaches the thermal 
radiation of the dS space. In the case of a contracting 
AdS universe we showed that the frequency of fluctu- 
ations produced at a certain moment grow in time and 
the temperature of radiation diverges as the big crunch is 
approaching, which is in agreement with earlier studies. 
We also explicitly demonstrated that when a mode of a 
certain frequency is created in dS space it does not auto- 
matically follow a thermal distribution with the temper- 
ature proportional to the value of the Hubble parameter. 
Modes whose frequencies are much larger or much smaller 
than the Hubble parameter do not follow a thermal distri- 
bution, while those of the order of the Hubble parameter 
do. Very high frequency modes are not thermal since 
they were produced from a non-thermal initial vacuum 
state and it takes some time for them to stretch and thcr- 
malize. Very low frequency modes which are larger than 
a horizon are frozen, and therefore non-thermal. As the 
high frequency modes get stretched by expansion, and 
their frequency becomes comparable to the Hubble pa- 
rameter, they become thermal. In AdS case, towards the 
Big Crunch, though the temperature nominally diverges, 
fluctuations of the spectrum grow and the temperature 
becomes less and less defined. 

The functional Schrodinger formalism we used here 
is different from the standard techniques used so far in 
similar studies. In Q, some objections to the standard 
quantum tunneling treatment (see for example [^-[IH) 
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are brought up, while in [8j it was suggested that the 
dS space is stable against the Hawking radiation. Our 
independent approach supports the standard picture of 
thermal radiation in dS space. 



z=0 



Therefore 



(A6) 
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in = M (i-^) n/2 n n (o). 



Since 



n n (o) 



( _ ir / 2 V2^^1, 



n=cvcn 



(A7) 



(A8) 



Appendix A: Number of particles radiated as a 
function of time 

We use the simple harmonic oscillator basis states but 
at a frequency ui to keep track of the different w's in the 
calculation. To evaluate the occupation numbers at time 
t > tf, we need only set Q = w(t/). So 



V2"n! 



-Hn{VmQb) (Al) 



where % n are the Hermite polynomials. Then Eq. ([22 
together with Eq. (TToT) gives 



1— 2 



1/4 e ia 



V2 n n\ 

1/4 e ia 



where 



P = 1 



V2"n! 



LJ \ p 



(A2) 



(A3) 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 



J(z) 



<%e- p t / 2 e~* + 2 ^ 



,-* 2 (l-2/P) 



(A4) 



and H n (0) = for odd n, we find the coefficients c„ for 
even values of n, 



(-l)" /2 e to [2f 2\ n/2 (n-1)!! 



For odd n, c„ = 0. 

We next find the number of particles produced. Let 



2 

1 -P 



(A10) 



Then 



N(t,u>) = 



n\c„ 



n—even 

2 d >r-^ 

V^\P\ X dx ^ 

v 1 1 ^ n—even 

2 d 1 

X 



(n-1)!! 



-A 



(1-x 2 ) 3 / 2 ' 

Inserting the expressions for \ an d P, leads to 



(All) 



N(t,u) 



ujp 

7f 



ujp 



(A12) 



Since 



(A5) 



In summary, we have found the occupation number of 
modes as a function of p which is a function of time as 
given by the non- linear differential equation Eq. (|17p . 
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